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Abstract 

We apply a recently proposed approximation method to the evaluation of non-Gaussian 
integral and anharmonic oscillator. The method makes use of the truncated perturbation 
series by recasting it via the modified Laplace integral representation. The modification of the 
Laplace transformation is such that the upper limit of integration is cut off and an extra term 
is added for the compensation. For the non-Gaussian integral, we find that the perturbation 
series can give accurate result and the obtained approximation converges to the exact result in 
the iV — >■ cx) limit (A^ denotes the order of perturbation expansion). In the case of anharmonic 
oscillator, we show that several order result yields good approximation of the ground state 
energy over the entire parameter space. The large order aspect is also investigated for the 
anharmonic oscillator. 



1 Introduction 

Anharmonic oscillator is a system that is well understood in both perturbative and non- 
perturbative aspects^'^. Still, the system plays an important role because it provides a good 
laboratory for the examination of any calculational scheme newly proposed. The Lagrangian 
of anharmonic oscillator is given by 

and the non-linearity is governed by the coupling constant A. The perturbation expansion is 
given in powers of X/nv'. For instance the ground state energy reads, 

E = mY.Aj^)\ (2) 



n=\) 



where the coefficients are found as^ 



1 . 3 , 21 

°"2' ^^4' '""T ''^''' ^^ 

It is known that the coefficient An grows as An ~ w ^(— 3)"+-'^r(n + 1/2) for large n (second 
paper in ref.l) and thus the series (2) diverges for any small A. Only for sufficiently small 
coupling constant, \/m? < 0.1, the series becomes numerically useful by the appropriate 
truncation. 

The linear S expansion^ is one of the framework which leads us to go beyond the weak cou- 
pling regime. For example, recent studies ^'^ succeeded to approximate the ground state energy 
in the pure anharmonic case {m = 0) by using the perturbative series (2). The strong coupling 
expansion of the ground state energy is also given in the literature^ (For the convergence issue, 
see ref.6). 

Recently, a new perturbative scheme was proposed in ref.7. The proposed method utilizes 
the information contained in the perturbation series by making Heaviside transformation^ with 
respect to some parameter (mass in ref.7). In the present paper we extend the method and 
explore how it can be used to obtain non-perturbative result in the two examples, a non- 
Gaussian integral and the anharmonic oscillator. We try to construct approximants for the 



integral and the ground state energy. For example, the approximant of the ground state energy 
will be given at A^-th perturbative order for A = 1 as 

_^.,.^ A„a;^"/^-V^ , ^ 7(3n/2 + l/2,mV) A^ 
^approxKru) e A.^p(3^/2 + l/2)^"';^o r(3n/2 + 1/2) m^-' ^^> 

where 7(p, z) denotes the incomplete Gamma function defined by 

7(p,z)= rdte-H^-\ (5) 

Jo 

and X* is the "cut-off" parameter to be determined order by order in some manner. By our 

approach one can approximate the ground state energy well over the entire region of the mass 

m, including the strong coupling limit and the crossover region from the weak to the strong 

coupling regime. We also note that application to field theories is straightforward. There 

exists no renormalization problem. 

There is a pioneering work by Graffi et al.^, where Borel summation method was used to 
compute the energy levels of the anharmonic oscillator. Our approach is similar to theirs. 
However, our approach is different from it in the respects that (i) the variable with respect 
to which the transformation is carried out is not the coupling constant but the mass square f\ 
and (ii) the Laplace integration is cut off and the extra term is added. We thus change the 
representation form itself, while ref.8 modified the integrand according to Fade construction 
to make up the Borel sum. In our case the integrand is just the truncated one. 

This paper is organized as follows: In the next section, we first review and then extend 
the method of ref.6. The extension results in modifying the ordinary Laplace integral rep- 
resentation as to fit the perturbative approach and enables us to deal with the effect of the 
explicit mass for both large and small m'^. In section 3 we perform a simple model calculation 
of non-Gaussian integral by using the method presented in the previous section. We show that 
the truncated perturbative series re-constructed by the method gives good approximation of 

the exact integral and converges to the exact answer in the N —>■ oo limit (A^ denotes the order 

tThc mass square, m^, corresponds to A^/"^ for the coupling constant. In terms of A^/'^/m^, the original 
perturbation series is not a power series. This is not a negligible difference in our scheme. The issue will be 
discussed in the last section. 



of expansion). In section 4 we turn to the anharmonic oscillator. Via perturbation series, the 
ground state energy is approximately calculated for various X/m? to several higher orders. We 
will show that, already at 5-th order, the error of our approximation is less than 1 percent 
for any \/m^ . In section 5 we address to large order aspects of our approach by proceeding 
to 249-th order. Some discussion and summary of the present work is given in section 5. In 
appendix we state the relation between our method and the linear 5 expansion in a particular 
limit. 

2 Modified Laplace transformation 

A part of the content of this section is same as the corresponding part in ref. 6. However, to 
make the presentation self-contained, we allow some overlaps with the work. 

For a given physical function /((t), we consider the Heaviside transform^ given by the 
Bromwich integral, 

rp+ico du \ 

Jp-ioD Zm a 

Here the parameter p represents the location of the vertical contour. Although rn^ corresponds 
to a in the cases we work with, there would be other choices in general. The contour of 
integration should be placed on the right of all the possible poles and the cut of f{a)/a. Then, 
if X < 0, the contour may be closed into the right half circle and f{x) is found to vanish. From 
f{x) we have /(o") via the Laplace integral of the second kind, 

dxexp{—ax)f{x). (7) 

-oo 

Since f{x) = when a; < 0, the integration range reduces to [0, oo). However it is convenient 
to keep the range as (— oo, +oo) to handle partial integration easily. 

Now suppose that one is interested in the value of / at cr = 0, but the perturbative expansion 
of /(o") does not allow one to let a arbitrary small. As in the Fourier transformation, small a 
behavior of /(o") is connected with the large x behavior of f{x). More precisely we find that 

hm /(a)= lim/(x), (8) 

o" — '■-f-U a:;— >oo 



where the existence of both hmits are assumed. Our approximation procedure is based on (8) 
and it goes as follows. Let /Ar(o') denotes the perturbative expansion of f{a) to A^-th order. 
Then the corresponding Heaviside function is given by 

fN{x) = / TT-^^fNi^)- (9) 

In cases of our interest, we can not take the naive limits, o" ^ or a; — > oo, in the both functions. 
Then there would be two routes to approximate /(O). Naive one is to approximate it by fixing 
a as small as possible in some manner. The other is, relying upon (8), to approximate /(O) by 
/jv(x) where x should be fixed at some large value, x* . Here note that f{x) often has larger 
convergence radius than /(a). For example, if f\a) = I]^o'^n/o"", then f{x) = J2'i^=o'^nX"'/n\ 
{x > 0). When f{x) has the convergence radius larger than that of /(cr), it would be convenient 
to deal with /^r rather than /jy. This is because we can probe the large x behavior of / by 
/at so that we have more chance to know the accurate value of /(oo) and /(O) accordingly. 
Therefore, we choose to approximate /(O) by /Ar(a;*). The explicit way of fixing x* is discussed 
in the next section. Here we just mention that, from the estimation of the upper bound of 
reliable perturbative region, x* will be fixed by the stationarity condition, 

dfN{x) 



0. (10) 



dx 

If there are several solutions, we should input the largest x* into f^. This is obvious because 
the value of / at a; = oo is what we are looking for. 

The above approach can be extended to the approximation of the function itself over the 
entire region of a. Let us start the discussion by showing how we can approximate the small 
a expansion of /(a), 

/(^) = /(0) + /«(0)a + /(2)(0)^ + ---. (11) 

As well as /(O), we can approximate f^''\0) as the following manner: From the formulas, 

df{a) df\x) 

o — -^ —X 



do dx 

-/(^) - r dyfiy), (12) 



a 



where the rightarrow represents the Heaviside transformation, we have 

f'\^) - £^ dy{-yf-M. 4.^ a,{x). (13) 

Assuming the expansion (11), the above two functions agree with each other at o" = and 
X = oo. Since what we have at hand is the perturbative one, /)/ (o"), the derivative 0:^(3;) is 
also truncated at order N. Hence, as in the previous case, we approximate f^'^'{0) by ak{x) 
by fixing the upper hmit of integration x as large as possible within the perturbative region. 
Replacing f{y) by fN{y), we then choose the input x, say xl, according to the same logic as 
that for X*. That is, x^. is determined by the stationarity condition. 



dxl dx% J -00 dy dxl 

We find that for any k 

xl = X* (15) 

where x* is a solution of (10). Substituting x* into a'^''\x) we can construct the approximate 
Taylor expansion, 

where aoi^x*) = fN{x*). 

Now it is an easy task to obtain expression which can be used for entire a region: We 
observe that the right hand side of (16) is easily summed to J^^ dxe^'^^df /dx. Then, by 
integrating by parts using f{x) = for a; < 0, it is written as 

e~"^*fN{x*) + a r dxe-^^^fNix)"^' Ma,x*). (17) 



The left hand side of (17) defines the modification of Laplace transformation. Note that, the 
naive x* — » 00 limit recovers the ordinary Laplace transform and gives /Ar(o"). We point out, 
however, that for small m the dominant contribution comes from the first term. Actually the 
second term, the cut-off Laplace integral, gives zero in the m ^ limit. Thus the first term is 
the crucial ingredient in our approach. To summarize, our approach results in approximating 
f{a)hj f^{a,x*) de&nedhy (17). 



It would be interesting to see how the approximant, fN{<^, x*), is different from the ordinary 
perturbative series. Let us rewrite the integral in (17) as 

a / dxe"'''fN{x) = /^{a) - a rfxe-"^/jv(x). (18) 



Then we have 



where 



/7v(a,x*) = fNicr) + f'^''''{a,x*), 



/r"(a, X*) = e-'^^ hix*) - a / dxe-^^fN{x). 

J X* 

If /Af(a) is given as /at = E^=o ^n/o-", then f^^"' is given by 

rcorr _ ^-(tx* \^ "'"•^ V^ 

In — ^ Z^ ^7j / . 

where 



^ a„a;*" A a„r(n + 1, ax*l 



n=0 



n! 



n=0 



n!(T' 



r(z,p) 

Using the asymptotic expansion of V{z,p) 
we find 



-t/2-l 



rfte"*t 



oo 1 

l + E-(^-l)(^-2)---(^-A:) 



I 



corr 

N 



_^^.x' g 1 d^fN{x*) 



k=l 



a^ dx 



*k 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



From its structure, (24) gives the large a expansion of f^^^. We find that under the stationarity 
condition (10) the first term in (24) vanishes. This shows that, for large a, the condition 
minimizes the deviation of the approximant from the ordinary perturbative result. 

3 A simple model calculation: non-Gaussian integral 

It would be worthwhile carrying out the calculation for a solvable model to gain some concrete 
feeling on the method to be used. We take up here a non-Gaussian integral. 



Z{m, A) 



2„2_5,„4 



-m'^q'^ — Xg^ 



(25) 



The perturbation series is given by expanding the integral in terms of A. The result reads 



oo 



n=0 "" '" 

It is easy to see that the series diverges for any \jvr&. Only when X/nv" < 0.1, the series can 
be put into numerical use by the appropriate truncation. Nevertheless, we will show that we 
can obtain the strong coupling {X/m^ ^ 1) expansion from the divergent weak coupling series, 
and even an approximant which is effective over the entire coupling regime. 

First step is to obtain the integrand of Laplace representation. This is done by calculating 
the integration (6) over m? (i,e., a = rn^). Then the result is given by 

°o ( — \\n ^2n+l/2 oo /_-|\n / /T \2n+l/2 

^ ' ^ to n\ (2n+l/2) ^ ' to n\ (2n + l/2) ^ '' ^ ' 

where d{x) denotes the ordinary step function {9{x) = 1 for a; > and for a; < 0). Note that 
Z converges for any Xx'^. For the sake of notational simplicity we set A = 1 hereafter. 
Let us introduce the perturbatively truncated Z by 

N (_'\\n 2n+l/2 

Z„(.,A=l) = Z.,W=gtfp-^^9W. (28) 

Now, we turn to the approximation of Z{m, A = 1) = Z{m) by using (28). We first discuss how 
to choose the input x* . The break down of perturbative series would generally appear as the 
rapid rise or fall of the Heaviside function due to the domination of the highest term in (28). 
Actually, the graphs of the function Zn{x) to the first several orders show it is the case (see 
Fig.l). Then, since Zn{x) is an alternative series, the rise and fall occurs alternatively. For 
example, for odd A^, Z]\f{x) temporally increases with x in the reliable perturbative region, but 
after that (eventually?) it falls down to — oo. We note that Z^r behaves temporarily flat just 
before the break down as shown in Fig.l. In this case the limit of perturbation would emerge 
typically as the stationary point. This is the reason we employ the stationarity condition, 

dx to "! „.o "! (2n+l/2) ^ '' ^ ' 

to fix the perturbative limit, x* {S{x) denotes the Dirac delta function). Dropping 6{x) and 



6{x) which are irrelevant, we have 

N t 1\n 

n\ 



oo 






ra=0 

The solution of (30) exists for odd A^ and depends on A^. By substituting x* into Zn{x) we 
have the approximant of Z{0). The obtained result is satisfactory as shown in Table 1. Note 
that X* increases with A^, which is a desirable result. 

We can also approximate the small mass (strong coupling) expansion of Z{m): First note 
that the coefficients of the k-th power of m^, ak{x*), is given by 

" ' J-^ ^ ' dx * ' ir„ n!(A + 2n + 1/2) ■ * ' 

Thus, at 15-th order for example, we have 

= 1.811655 - 0.609988m2 + 0.223363m^ - 0.074001m^ + 0.022042m^ + ■ ■ (32) 

The exact result reads from (25) that Z{rn) = |l]^or(^/2 + l/4)(— 'm^)"'/n! and is given 
numerically as 

Z{m) = 1.812805 - 0.612708^^ + 0.226601m^ - 0.076589771*^ + 0.023604m® + ■ ■ ■ . (33) 

We see the good agreement of (32) and (33) up to several orders. Thus, we reach sufficient 
accuracy of strong coupling expansion from only the information of the truncated weak coupling 
series. 

Next we examine the approximation of Z{m) for various m. The approximant is given from 
(28) and (17) as, 

7 (m r*) - P-^^*^ f i-^ri^rm^ f (-1)" 7(2n+l/2,mV^) 
"^^ '^" h n\ (2n + l/2)+„^,n!(2n + l/2) ^^^^ ' ^^'^ 

Table 2 shows the value of the approximant at A^ = 15 together with the exact results. We 

find that the approximant gives good values for the sample of m. 

Before closing this section we show that the approximant Z]^{x*) converges to Z{oo) in the 

N —^ oo limit. First we show the convergence of lim2;^oo Z{x)^ where Z{x) denotes the exact 



transformed function, by calculating Bromwich integral exactly. From (6) and (25) Z{x) is 

given as 

Z{x) = 2 f^dqe-^' = 1-^(1/ A, x'), (35) 

Jo 2 

where we used that e'"^ '^ transforms to 6{x — q^). Thus it is apparent that lim2:_+oo Z{x) = 
r(l/4)/2, which agrees with Z{0). 

For the proof we need to know how x* behaves for large N. The relation is found as follows: 
We note that the condition (30) is viewed as the truncation of the equation, x~^/^ exp(— a;^) = 0. 
Since the series expansion of exp(— x) has infinitely large convergence radius, the obtained 
solution tends to +oo as A^ ^ oo. More precisely, by assuming the form, x*'^ ~ aN^ (a, b: 
constant), we find the following scaling at large A^, 



*2 
X ~ 



1 



-N\ (36) 

Now, let us define the reminder, R^, by 

Since Zoo{x), the perturbative series to all orders, apparently converges to exact Z{oo) because 
of the infinite convergence radius, it is sufficient to show that 

lim Rn{x*) = 0. (38) 

Af^oo 

This is easily verified: Using the Stiring's formula, we obtain 

'^ „n/ *2\n+l/4 



and from (36) we then find 



\RN{xn\<^^lf^N-'/'Q''^0 (iV^oo), (40) 

Vo7r(l — e/3) ^3^ 



which proves (38). 
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4 Approximation of the ground state energy of the anharmonic 
oscillator 

We turn to discuss the anharmonic oscillator from this section. In complicated systems it is 
generally hard to proceed to arbitrary higher orders. Hence it is practically important to study 
whether an employed approximation scheme works at low orders. 

In this section we use perturbative series up to 9-th order and show that the several low 
order result can yield good approximation of the ground state energy via modified Laplace 
representation. 

Let A = 1 for notational simplicity. Our first task is to obtain the Heaviside function of the 
perturbative ground state energy, Eiy{m), 

EnM = g ^^2)3n/2-l/2 - (^1) 

Heaviside transform of rn^'^"+^ with respect to m^ gives 

3n/2-l/2 

m-3"+i ^ ——, r^dix). (42) 

r(3n/2 + l/2) ^ ^ ^ ^ 

Then from (41) and (42) we have the Heaviside function, Efyf{x), 

N A 3n/2-l/2 

^"'-'^g. ^(3;/2.l/2) ''W■ («' 

For notational simplicity we omit the step function in what follows. 

By suitable replacement of variable, the function E agrees with the function appeared in 
ref. 5 in which the 6 expansion method was applied to the anharmonic oscillator (see eq.(7) of 
ref. 5). Authors of ref. 5 considered the case where the coupling, SQ"^, grows to large orders {6 
represents a fictitious parameter which is to be set 1 at the end of calculation and Q^ denotes 
the mass auxiliary introduced to divide the given Lagrangian un-conventionally) . Then the 
authors found that their approximant converges to the form (43) in that large order limit 

0. This implies that there may be some connection between the linear 6 expansion and our 

tin S expansion scheme, this Uniit corresponds to infinite perturbative order. Then we note that in this 
hmit the parameter fl is no longer arbitrary because it should be of the same order of magnitude with the 
perturbative order. We also point out that this limit was considered in ref.9 from different motivation. 
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method. Actually we show in appendix that it is the case; We will show that, in a particular 
limit, the Heaviside transform is induced in a suitable interpretation. 

Now, let us discuss the energy approximation using (43). As in the previous section we use 
the following condition, 

^^0, (44, 

to fix the input x*. Up to A^ = 9, while there is no solution for even A^, we find just one 
solution for each odd A^. By substituting x* into En{x), we obtain the approximation of E{0) 
which is the energy in the pure anharmonic case or in other words in the strong coupling limit. 
We also evaluate the succeeding Taylor coefficients of first five terms. For the purpose it may 
be convenient to use the integrated form of ak, 

^ ^(0) ~ «, = (-1) E ^^3^^2 + l/2)ik + 3n/2 - 1/2) (^ = 1' 2, 3, ■ ■ ■), (45) 

which comes from substitution of (43) into the general formula (13). Then, we have the strong 
coupling expansion, 

E{m) ~ ao + ^m^ + |^m^ + ■ ■ • , (46) 

where ao = En{x*). We have done numerical calculation with Mathematica. As shown in 
Table 3, the obtained results agree well with the recent result reported by Kleinert^. Next we 
check the approximation for various m. For the purpose we substitute (43) into (17) and use 
the definition of incomplete Gamma function (5) to obtain the approximant, EN{m,x*). The 
result reads at the order A^ as 

£,„(,„.., 2_.r(3„/2 + l/2)+'"„% r(3n/2+l/2) m»»- *""' 

We find that the numerical calculation shows the good results for all m |. From these results, 
we find that the accuracy of calculated energy is quite satisfactory. Already at the 1-st order, 
the approximation gives error only within 11 percent for all m. And at 5-th order, the error is 

less than 1 percent. These results are depicted in Fig. 2. In particular we note that E]\;{m, x*) 

§The reference value is generated according to the method of ref.ll. We thank H. Suzuki for the correspon- 
dence of the reference. 
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improves the approximation of E in the crossover region of weak and strong couphng regimes, 
0.1 < 1/m'^ < 1. Note that in this intermediate couphng region the ordinary perturbation series 
can not be used because the series never be close to the exact value under any truncation. 

It would be interesting to see the explicit difference of our approximant from the ordinary 
perturbative series. Repeating the steps from (18) to (24), we find 

Er,{m,x*) = EN{m) + E^°''^ (m, x* ) , (48) 

where 

E^ [m,x)-e E^[x) m 2^^ r(3n/2 + 1/2) m^" " ^ ^ (m^)^' ^^^^ 

and 

^0 r(3n/2 + l/2) /i/ ^ ^ -^^ 9x« ^ ^ 

The coefficients slowly changes order by order except for 61 jv which satisfies from (44), 

bi,N = 0. (51) 

We have done numerical computation of bi^N at A^ = 1, 3, 5, 7, 9. To several i, we find that the 

size of bi^M decreases as N increases. This is an expected tendency. In fact, we can show that 

bi^isr should converges to zero in the N —* 00 limit. The basic relation to be noted is that 

dfjx) 
af{a) -. ^-. (52) 



From (52), it is easy to see that 



and for our case. 



^'/(.)^^, (53, 



m'^Eim) ^ ^13^, (54) 

^ ^ dx' ^ ^ 

There should be the agreement condition (8) between these functions and therefore, noting 
that lim.fn2-,o'nT''^^E{m) = 0, the coefficients bi^^ should tend to zero as A^ increases to infinity 
if the approximation procedure is working well. This issue will be studied further in the next 
section. 
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5 Higher order behavior 

The important information in our approach is contained in the Heaviside function E{x). We 
therefore investigate the properties of E{x) relevant to our analysis by extending the pertur- 
bative order up to 249-th. 

An important issue in our approach is whether lim^_^oo E{x) exists or not. Although we 
do not have rigorous proof, we see convincing answer by figuring out En{x) to large A^. We 
have generated perturbative coefficients An up to 249 terms with the help of Mathematica and 
plotted the graph of -^249(2;) as shown in Fig.3. We note that plateau starts around x ~ 1 and 
abruptly grows up around x ~ 3.2, which shows the break down of perturbation expansion. 
Taking closer look, we find that the function E24g{x) weakly oscillates at the plateau region. 
The amplitude of the oscillation is very tiny indeed; The difference between the first extremum 
and the next is just 0.0000103 ■ ■ ■ which should be compared with the first extremum value, 
0.667975902279 ■■ ■. The difference between the second and third is about 0.00000001. Thus 
the amplitude decreases as the function oscillates to larger x. The values of three stationary 
points of -E249 are given as 

E249 = 0.667975902279 ••• , a; = 1.139689002700 
E249 = 0.667986268727- ■■, a; = 2.069065340532 
E249 = 0.667986259143 ■■•, a; = 2.987637042160. (55) 

These values shows how the En' at the plateau is close to the value E{0) which is known to 
bei2 

E{0) = 0.667986259155777108270962 ■ ■ ■ . (56) 

Thus the behavior of ^249 (2^) for a: < 3, where the function is reliable, strongly suggests that 
lim.x-,00 E{x) would exist and consequently agree with E{0). 

The first, second and third stationary points begins to appear from 28,101 and 246 orders, 
respectively. Hence, several solutions exist for some orders higher than 28-th. This phe- 
nomenon was also observed in the 6 expansion framework^. There it was found that the value 
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of interest at largest l/Vt was most accurate. But within the framework there is no a priori 
reason why one should take the largest l/Vt. On the other hand, it is obvious in our approach 
that one should focus on the largest x* as we mentioned in section 2. The behavior of our 
approximants En{x*) as the order increases is as follows. For x smaller than the largest x* , 
En{x) is a good approximation of the exact function and the departure starts around x ^ x*. 
Therefore the largest stationary point sits in the vicinity of the exact function, and slides to 
larger x direction along the curve of the exact function as N increases. Hence the convergence 
issue of the approximants, {E]\f{x*)\N = 1,2, ■ ■ ■}, is tightly connected with the convergence 
of linia-^oo E{x) and how the stationary solution, x*, grows with the order. Since, by the def- 
inition, the largest x* is located at the upper limit of reliable region, {En{x*)\N = 1,2, ■ ■ ■} 
would converge to E{oo). It is also clear that En{x*) oscillates as A^ increases by following the 
function E{x). This oscillatory property of the approximant was observed (but not clarified) 
by Kleinert^. 

As is obvious from the above discussion, the largest stationary point around the order 
A^ = 246 (at which the third stationary point of E{x) is settled) is approximately given by 
the third stationary point shown in the last of (55). Therefore using that value of x*, we can 
see how accurate Ei\i{m,x*) is for various m^. The result of computer calculation is shown in 
Table 4 and shows that the obtained values are quite accurate for all m^ . 

Finally let us comment on the behavior of fej^jv, the coefficients of E'^^^ at large m, for large 
A^. We have calculated them at A^ = 28, 101, 249. We find that the size of bi^N decreases to zero 
as the order A^ increases. For example, the results for i = 1 to 7 at A^ = 249 are respectively 
given as follows; 

6i,249 = OE-22, 62,249 = 8.259931E-10, 63,249 = -8.27746E-9, 64,249 = 5.094257E-7, 
65,249 = 4.804239E-5, 60,249 = 6.054357E-3, 67,249 = 7.451039E-1, (57) 

Thus, for lower i, the result almost agrees with the requirement that lim.x^ood^E{x)/dx^ = 0. 
For larger i, however, the coefficient 6j 249 grows rapidly. This represents that as x exceeds x* 
E2A9{x) abruptly grows as shown in Fig.3. 
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6 Discussion and conclusion 

Critical things leading our scheme to the success are that the convergence radius, p, is infinite 
for the Heaviside functions f{x) and that they quickly approach to the value f{oo) at finite 
X. The later fact is confirmed numerically for the anharmonic oscillator from the remarkable 
closeness of En{x) at a; G [1,3] to E{oo). For the non-Gaussian integral it is analytically 
clarified as follows: By using 7(1/4, x^) = r(l/4) — r(l/4,a;^) and the expansion (23) we 
obtain from (35) that 

1 



Z{x, 2 



r(l/4) - a;-3/2 exp(-a;2) 1 + E ^"' IT ( V4 

^ i=i i=i 



(5^ 



The approach to /(oo) is speedy because of exponential damping at large x. It is interesting 
to note that Z{x) has no power-like term, a;~", in large x expansion (58). 

Actually these features are related with the choice of integration variable generally denoted 
as (7. We discuss on this issue by our two examples. 

First consider the non-Gaussian integral case and let a = m^. We numerically study on 
how the function Z varies according to the power f3. In calculating the transform, the following 
result is convenient to use: 



a 



e 



X 



-« 



r(i-0' <''' 

where we have dropped possible 6 and 6 functions. From (59) it is easy to find that, with 
respect to m^, the transformed series is given by 

Z^(x) = Vf-ir r(2n+l/2) 4n//3+l//3 (QQ^ 

Then we find that Zoo{x) is a divergent series for (3 > 4 and becomes a convergent one for 
(3 = 4 where the convergence radius is 1/4. When (3 < 4, the series becomes a convergent one 
for any x. We have plotted in Fig.4 the series for j3 G [1,3] at A^ = 100. Next problem is 
whether the limit lim.x-*oD Zoo{x) exists or not. From the numerical calculation at A^ = 200, 
we find that the result changes around j3 ~ 1.32: For j3 G [1, 1.32], ^00(2^) would not converge 
in the a; ^ 00 limit. The convergence of lim^i^oo Zoo{x) seems to be realized for (3 G [1.32,4]. 
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We do not know why around [3 = 1.32 the convergence property changes. Now, we observe the 
following features: When (3 is small the truncated series becomes effective to large x but Z{x) 
does not converge at a; = oo. As (3 increases beyond ~ 1.32 but still within /3 < 4, the effective 
range of truncated series decreases but the saturation to finite Z{oo) is fast. Furthermore the 
oscillation amplitude damps as (3 decreases. From these points, it is subtle that, at a fixed 
order, what value oi [3 between ~ 1.32 and 4 gives the best approximation. We have calculated 
Zn{x*) at N = 99 or 100 for (3 = 1.5, 1.7, 1.9, 2.0 and 2.1. The result reads 

p = 1.5 1.8169575 {x = 15.037, N = 100) 

p = 1.7 1.81926459 {x = 10.336, A^ = 100) 

p = 1.9 1.81719914639 {x = 7.1806, A^ = 100) 

p = 2.0 1.812804954110934 {x = 5.34, A^ = 99) 

p = 2.l 1.8051679 (a; = 5.043, A^ = 99). (61) 

Thus we find that the best approximation is realized at /? = 2.0. The agreement to the exact 
value, 1.812804954110954 ■■ ■ is remarkable for that case (The first 14 decimals are in the 
agreement). 

We note that the ordinary Borel summation method employees A as the transformation 
variable. This means in our context that a = m'^. Then the convergence radius of pertur- 
bative Z{x) becomes finite (p = 1/4) and the choice is not suited for our modified Laplace 
representation approach where the integrand is truncated and the upper integration limit is 
cut off. 

Now we turn to the anharmonic oscillator. If one performs Heaviside transformation with 
respect to m^, one has 

N A 

En{x) = E fTFTTT^^TTTT-TT^^'""'^^'- (62) 

^0 T{3n/P - 1/P + 1) 

The n-th coefficient behaves for large n as 

Ti3n/P-1/P + 1) ^ ^ \ej ^ ^ ^ 
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Hence the series is divergent for f3 > 3 and p is infinite for f3 < 3. For (3 = 3 the series is 
a convergent one, but p = 1/3. We have figured out -^249(2;) by varying (3 from 1 to 3. The 
resuh is shown in Fig.5. The convergence property of Eoo{x) in the a; ^ 00 hniit seems to 
change around (3 ~ 1.15. As in the previous case, we examined what value around (3 = 2 gives 
the best approximation of E{0) at N = 248 or 249. For (3 = 1.7, 1.8, 1.9, 2.0, 2.1 we obtained 

p = 1.7 0.6655 {x = 6.53 ~ 6.57, N = 248) 

/9=1.8 0.665514724 (x = 5.1257 ~ 5.1260, A^ = 248) 

p = 1.9 0.666234824085 (x = 4.0185, A^ = 248) 

p = 2.0 0.667986259143255939 {x = 2.98685, A^ = 248) 

p = 2.1 0.67107970759 {x = 2.5012, A^ = 249). (64) 

From the above sample, we find that the choice /3 = 2 gives the particular good approximation. 
The ordinary Borel choice of variable, a = m^, also does not work well in our approach. This 
is because in that choice the convergence radius of perturbative E{x) becomes finite (p = 1/3). 
It is interesting to note that why the choice, P = 2, gives the best approximation in both 
cases. Although we have not resolved it yet, we suspect that the resolution would have related 
to the fact that Z{m) and E{m) has expansion in the square of the mass, Cq + Cim^ + 02(^72^)^ + 
• • -. Any finite sum, J2^=o Ci{'m'^y , vanishes, if the integration variable is chosen as cr = m^ 
because from (59), 

a'^0 (i = l,2,---). (65) 

Hence, there is no powers of 1/x in the corresponding Heaviside function at large x. Note that 
when the power ^ is fractional the transform of cr^ survives to give x^^. Thus, as well as Z{x) 
which is analytically solved as (58), E{x) would be suppressed exponentially only at /9 = 2 
and the smallness of the difference \E{x) — E{oo)\ would be thus explained. 

To conclude this paper, we have demonstrated how truncated perturbation series can be 
utilized to approximate the non-perturbative quantities via modified Laplace representation. 
The non-Gaussian integral is precisely calculated from the truncated perturbation series. And 
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the approximant for the strong couphng hmit is found to converge in the N ^ oo hmit. 
The ground state energy of the anharmonic oscillator is also approximated successfully over 
the entire parameter space. Large order calculation strongly suggests that our approach works 
well to all orders. In both examples, we find that the perturbative knowledge serves us sufficient 
information for recovering the small mass or strong coupling expansion. 

In field theories, however, perturbation series has more complex structure in general. For 
example, the perturbation series for the dressed mass, vacuum energy density or condensates 
involves the mass logarithm, logm. When the mass-logs exist, it is no longer expected that the 
convergence radius is enlarged by Heaviside transformation with respect to the mass, although 
some qualitative improvement is found in the literature (first paper in ref.6). This is the 
subject of our future investigation. 



19 



Appendix: Linear 6 expansion in a particular limit and the Heaviside 
transformation 

In this Appendix we show that, in a particular hmit, the hnear S expansion leads to the 
Heaviside transformation. 

In linear 6 expansion, the Lagrangian of anharmonic oscillator is written by introducing 
auxiliary or variational mass Q as, 

and the free and interaction parts are defined as 

Lfree = 7T ( 37 ) ~o^ 1 ' 



2\dtJ 2 
Lint = ^5n\'-\q\ (67) 

where 

S = l- myQ\ (68) 

Let us concentrate on the 171 = case hereafter. Then one should set 5 = 1, which represents 
the masslessness of the starting Lagrangian, at the end of the calculation. 

First we point out that perturbation expansion in Lint is generated from the ordinary 
perturbative result with mass Q by shifting Q"^ ^ Q'^{l — 6) and then expanding in powers of 6 
to relevant orders. Let the number of SQ"^ vertex and A coupling in a given Feynman diagram 
be no and n\ respectively. Then at order n in Lint, any contributing diagram should obey 
nn + nx < n. To proceed to large orders keeping the analytical manipulation straightforward, 
it is however convenient to modify the expansion scheme. The new expansion is defined such 
that, at order n, contributing diagram should obey n\< n and rin < n. Even a diagram has n 
vertices for A coupling and n mass insertions, it should be included to the n-th order expansion. 
In this expansion scheme, any contributing diagrams are obtained such that, given a Feynman 
diagram with no 5Vl? vertex, one should incorporate it by shifting 5 and then expanding in 5 
up to just the order of perturbative expansion. This procedure is formally carried out as we 
can see below. 
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Consider a given Feynman amplitude with no 6Q'^ vertex, /(fi^). First we shift Q^ 
r2^(l —6) and expand the result in 6. We then have 

na'-m, = f:fi^>(n^-)t^ = H„j_f; (z^(_4)'=;(sj=). (69) 

Setting 5 = 1 we have 



k\ N^^^^ k\ ^dn^ 



N (_o2\k Q , 

m = }ffiT^-^(^)m% (70) 



By using 



we find 



Thus we arrive at 



^^'^t^o k\ V91]2 



9 1 _ 1 9 1 
dpp pdp p^ ' 



(71) 



1^^:^[7^) =P^7^[l^) - = ^n{p). (72) 



fe=0 



k\ ^dp' N\ ^dp' p 



/(O) = hm P^(l]2)/(i]2)_ (73) 

This is, however, just a formal result and we address explicit example and ask what comes out. 
The operation of ^^{Vt^) on (f2^)^ gives 

When N is large enough, the right hand side approaches to (r2^/A^)^/r(l — ^) and the result 
of operation of ^^{Vt^) on E^i^) reads, 

^"'"''^-'»)-|: r(3W2;i/2) (^)'"'"' (^»^'- P^' 

Thus when N/Vi^ is changed to x in (75) the result agrees with (43). However note that, in 
accord with iV, we must let Vt large enough to stop the divergence of RHS of (75) (Note that at 
N = 249, the largest stationary point is given at N/Vt^ ~ 3). To give (75) some meaning, it is 
thus necessary to let A^ and Vt^ simultaneousely large with N/Vt^ kept finite. And accordingly, 
Vl? is no longer a free finite parameter. If we adopt this recipe, the equality of (73) breaks 
down and we should take the righthand side of (73) as a new function / of variable x = N/D? . 
We can expect the equality relation only in the limit, x — * oo, such that lim^^oo f{x) = /(O) 
(This limit corresponds to A^ ^ oo with ^2^ fixed). 
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Now we show that V^ induces Heaviside function in the hmit that iV, f2 ^ oo with Q'^/N 
fixed. Let us start with the Laplace transform, 

/oo 
dtexp{-nH)f{t). (76) 

-oo 

By operating Vn we have 

/■oo (n2\N+lj.N „ , , /.oo 

Vr^fin) = j Jt^—L^^^ eM-^'t)f{t) = j JtAN,a.{t)f{t). (77) 

We find that the function Ajv,Q2(t) approaches to the Dirac 5 function in the hmit, A^, ^2 ^ oo 

with the ratio N/Vl? fixed, 

A^,f,2(t)^5(t-iV/fi2). (78) 



Thus we find that 



Jim P^/(fi) = /(iV/fi^). (79) 

7V,C2— »oo 



This result states that T>n gives the kernel of Laplace transform in the limit N,VL ^ oo with 
the ratio N/Vt^ fixed. 

We remark that although we have shown that, in our modified expansion scheme, the 5 
expansion can induce the Heaviside transform in large orders, this may not lead that the 
conventional expansion in Lint also gives the same function / in the large order limit. This is 
because the difference of the included diagrams between the two expansions increases as the 
order increases. 
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Table Captions 
Table 1 Numerical computation of Zj^{x*) to 15-th orders. 

Table 2 Numerical computation of Zi^{m? ,x*) for various values of m^ . At this order x*"^ = 
5.0438870. 

Table 3 Low order approximation of first several strong coupling coefficients of the ground 
state energy. 

Table 4 Numerical result of E24g {m^ , x* ) for various rn^ . 

Figure Captions 

Fig. 1 The perturbative Heaviside functions Ziq{x) for A^ = 1,4, 7. 

Fig. 2 The ratio, E]\f{m, x*) / Eexacti'nT') is plotted at A^ = 1 and 5. 

Fig. 3 The function E24g{x) is shown. The sharp rise around x ~ 3.2 represents the break 
down of the perturbative truncation. 

Fig. 4 The function Ziqq{x) is shown for various choice of (3. The amplitude of oscillation of 
Ziqq{x) decreases as f3 decreases from 1.0 to 2.0. The effective range of the series is larger 
for smaller f3 but the saturation to Z{oo) is slower. When f3 is close to 2.0, the oscillation 
is weak and the saturation to Z{oo) is very fast. 

Fig. 5 The function -^249(2^) is shown for various j3. At /? = 1, En{x) ^ as a; ^ 0. When 
1 < /5 < 3, En{x) diverges in the x ^ limit. The oscillation amplitude of E2'i9{x) 
damps as f3 decreases. At f3 = 2.0, the oscillation is quite weak and the saturation to 
E{oo) is very fast. 
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N 


Zn{x*) 


x*2 


1 


1.6 


1 


3 


1.7313594 


1.5960716 


5 


1.7765256 


2.1806071 


7 


1.7955618 


2.7590027 


9 


1.8043006 


3.3335514 


11 


1.8085078 


3.9054517 


13 


1.8105959 


4.4754119 


15 


1.8116546 


5.0438870 


exact 


1.8128049 





Table 1 



m 



Zi^{m\x*) 



exact 



0.01 1.80557702921362946 

0.1 1.75281762908207768 

1 1.36831695165151724 

3 0.96173724333279584 

6 0.710038679143677110 

10 0.556465718382570615 

100 0.177232097497741759 



1.806700454307384679 
1.753725831772014832 
1.368426855735508774 
0.961738333157472108 
0.710038680405132855 
0.556465718382772753 
0.177232097497741761 



Table 2 



N 



E(x*] 



«i 



^2 



1 0.738558766382022 

3 0.686283726385561 

5 0.674564660427775 

7 0.670682699394559 

9 0.669175108154224 

exact 0.667986259155... 



0.121215344755496 
0.134984882799344 
0.139821127686036 
0.141853009994811 
0.142780937393257 
0.143668783380... 



-0.004420970641441 
-0.006437328047158 
-0.007465973798930 
-0.008005281363443 
-0.008293166786158 
-0.008627565680... 



N 



"3 



^4 



X^ 



1 0.0002176764415436 

3 0.0004280233911678 

5 0.0005773467460109 

7 0.0006739106832314 

9 0.0007338756445212 

exact 0.000818208905... 



-0.000010207418675 
-0.000027447789756 
-0.000044015663387 
-0.000057114935211 
-0.000066541377300 
-0.000082429217... 

Table 3 



0.328248340614232 
0.448360373271548 
0.549152913559036 
0.636621775859320 
0.713814879924458 
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m 



E2A^{m\x*] 



exact 



0.001 0.66812991929974827 

0.01 0.66942208503810206 

0.1 0.68226767187380087 

1 0.80377065123375873 

10 1.64938954183035211 

100 5.00747395574729234 
1000 15.8121382178529 

Table 4 



0.66812991931241042 
0.66942208505040309 
0.68226767188301217 
0.80377065123427376 
1.64938954183035211 
5.00747395574729234 
15.8121382178529 
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